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It is found that the zero field susceptibility \ of the Ising model on the Cayley tree exhibits unusually 
weak divergence at the critical point Tc- The susceptibility amplitude is found to diverge at Tc 
proportionally to the tree generation level n, while the behavior of % ls otherwise analytic in the 
vicinity of Tc, with the critical exponent 7 = 0. 
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It has been well established^ □ that the Ising model 
on the Cayley tree exhibits no spontaneous order (zero 
magnetization in zero field), while the susceptibility di- 
verges for temperatures lower than the critical value 
T c = 2fc B 1 J/ln[(VS + l)/(y/B- 1)], where k B is the 
Boltzmann constant, J is the nearest-neighbor interac- 
tion parameter, and B is the tree branching number (co- 
ordination number minus one). Above Tc, it was [found 
through several different approximation schemesa a that 
the susceptibility asymptotically behaves as 



/3{u+l) 2 
(l-Bu 2 



(1) 



where (3 = 1/ksT and u = tanh(/3J). Expanding the 
above expression around t= (T — Tc)/Tc one obtains 



attention to its singular behavior in the thermodynamic 
limit, at and below Tc- In particular, it turns out that 
the singular behavior of x a t Tc , as displayed by the ex- 
act formula, is quite different from that deduced from the 
appropcixaate asymptotical expression (Q) used in previous 
worksDtira. More precisely, we find that the terms of the 
order ~ t~ l exactly cancel out, and it is shown that the 
susceptibility displays what may be termed a "divergent 
coincident singularity", in Fisher's classifications, with 
critical exponent 7 = 0. 

For simplicity, we consider here only the Cayley tree 
with the branching number B — 2, the generalization to 
arbitrary B being straightforward. We further consider 
a single n-generation branch of a Cayley tree, composed 
of two (n — l)-generation branches connected to a single 
initial site, with the Hamiltonian 



Vb + i 



2jVb (Vb-i 



t- 1 + O (t°) 



(2) 



n 



-J SiSj — H Si 
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corresponding to the "classical" value of the critical 
exponent 7 = 1. The same value 7 = 1 was also 
foundBQ for the Bethe lattice (interior of the Cayley 
tree), with the (higher) critical temperature Tbc — 
2k B 1 J/\n[(B + l)/(B-l)}. 

On the other hand, in a recent workQ the present au- 
thors have established the exact closed form expression 
for the zero-field susceptibility x { n : T) of the Ising model 
on a Cayley tree of arbitraiju.|generation n. Our resultLl 
confirmed previous findingsaljlj concerning susceptibility 
high temperature behavior (T > Tc), but it also sheds 
new light on its singular behavior in the low temperature 
region T < Tc ■ In this work we perform a finite size scal- 
ing analysis of the zero field susceptibility, paying special 



where H is the external magnetic field, Si — ±1 is 
the spin at site i, (nn) denotes summation over the 
nearest-neighbor pairs, and J denotes the coupling con- 
stant as before. The n-generation branch consists of 
N n — 2 n+1 — 1 spins, the 0-generation branch being a 
single spin. The recursion relations for the partition 
function of any twexconsecutive generation branches were 
found by EggarterO to be 



7 ± - i^ 1 



(4) 



where x = exp(/3J), y = exp(/3H), and Z+ and Z~ de- 
note the branch partition functions restricted by fixing 
the initial spin (connecting the two (n — l)-generation 
branches) into the {+} and {— } position, respectively. 
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As the nonlinear coupled recursion relations (Q) can 
be iterated to yield a closed form expressiona only in the 
zero field case n several rather sophisticated approxima- 
tion schemesoQ have been devised to elaborate the field 
dependence of the partition function, and deduce (among 
other quantities of interest) the limiting behavior of the 
zero field susceptibility in the high temperature region, as 
given by equation ([!]). _ 

On the other hand, only recentlyQ we have arrived at 
the exact analytical expression for the zero field suscepti- 
bility, by considering the recursion relations for the field 
derivatives of the partition function, which can be iter- 
ated in the limit H — ► to yield corresponding closed 
form expressions. Zero field susceptibility of a Cayley 
tree of arbitrary generation n, as a function of tempera- 
ture, was thus foundQ to be 
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oo the first term (h n ) inside the 
square brackets on the right hand side of (||) recovers 
formula (Fy), and represents the dominant term for tem- 

(2) 

peratures T > Tc- The second term {h n ) is dominant 
in the low temperature region T < Tc, while the last 
two terms (h^h 1 and h^) become negligible in the ther- 
modynamic limit for all temperatures, and are therefore 
relevant only for finite size systems. 

Up to this point everything matches-the previous re- 
sults and conclusions of other authoraJa, therefore, it 
came as somewhat of a surprise to find that upon expan- 
sion in power series around t = the terms proportional 
to t^ 1 exactly cancel out, independent of the system gen- 
eration level (therefore, also in the thermodynamic limit), 
on both sides of the critical point. It will be shown in the 
remainder of this paper that this fact leads to quite differ- 
ent conclusions about critical behavior of Xi from those 
that follow from the approximate expression ([!]). 

Expanding the four individual sequential terms, fin , 
hn j hn^ and hlt\ inside the brackets on the right hand 
side of ([5]), in the vicinity of Tc, and retaining only the 
terms of the order it is found 
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where K c = J/k B T c = In (l + ^2) ~ 0.881374. It is 
seen that the terms proportional to t^ 1 exactly cancel 
out for arbitrary n, for both positive and negative t. Re- 
taining additional powers in t, for large n one obtains 



Xn = nK c 



2J 



1 - 1 n K c -= + 0{t 2 ) 
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which is valid on both sides of the critical point for 
\t\ « 1/n (since t is an independent small parameter, 
this condition can be fulfilled for arbitrarily large n). It 
is seen from (0) that in the thermodynamic limit Xn di- 
verges at Tq proportionally to tree generation level n, 
rather then as t^ 1 as follows from approximate formula 
(|l|), and is otherwise analytic in t. It therefore follows 
that the susceptibility critical exponent is equal to 7 = 0, 
while the amplitude diverges as n — ► 00, demonstrating 
behavior that may be termedu "divergent coincident sin- 
gularity" . It should be noted that a Cayley tree of gen- 
eration n = 77 has a number of spins corresponding to 
the Avogadro's number, while for n ~ 270 the number of 
spins corresponds to the estimated number of hadrons in 
the observable Universe. Therefore, the above suscepti- 
bility exhibits an extremely weak singularity. 

To illustrate the behavior of susceptibility in the vicin- 
ity of Tc in Fig. 1 we show the function Xn/n, as given by 
formula (^), for several system sizes n = 64, 128, 256, 512 
and 1024. All the shown curves Xn/ n intersect at 
k B T c /J = l/ln(l + V2) - 1.134593, exhibiting a fi- 
nite value In (l + V2) (l + \/2) 2 /2 - 2.568511. Note 
that these correction terms turn already imperceptible 
on the scale of the graph. 




FIG. 1. Zero field susceptibility as a function of temper- 
ature, calculated using formula (Bl), for several system sizes 
n = 32, 64, 126, 256 and 1024. Crossing of the curves /n at 
Tc demonstrates the unusually weak singular behavior. 

In the temperature region below Tc the second term 
(hn^) inside the brackets on the right hand side of (|5|) 
becomes dominant, and for large n susceptibility is given 
by 



Xn= P- 



u 2 (2u 2 ) 



n+l 



(2u 2 - l)(2u- l) 2 ' 
It is seen that the susceptibility diverges as (2u 2 



(8) 
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the divergence becoming stronger as temperature is low- 
ered. Scaling of susceptibility for different system sizes 
is obtained by taking the logarithm and then dividing 
by the generation level n. To demonstrate this scaling 
behavior, in Fig. 2 we show the function lnx„/ (n + 1) 
for several system sizes, together with the corresponding 
limiting function In (2m 2 ) . It is seen that the displayed 
finite size curves become hardly distinguishable from the 
limiting function for < T < Tc- 
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FIG. 2. Finite size scaling of the zero field susceptibil- 
ity in the temperature region < T < Tc- The shown 
curves are calculated using formula (||), for several system 
sizes n = 32, 64, 126, 256 and 1024. The dotted line repre- 
sents the limiting curve In [2tanh 2 (/3J)l. 



1 J. von Heimburg and H. Thomas, J. Phys. C 7 (1974) 3433. 

2 H. Matsuda, Prog. Theor. Phys. 51 (1974) 1053. 

3 E. Muller-Hartmann and J. Zittartz, Phys. Rev. Lett. 33 
(1974) 893. 

4 T. Morita and T. Horiguchi, J. Stat. Phys. 26 (1981) 665. 
R.J. Baxter, Exactly solvable models in Statistical Mechan- 
ics, Academic Press Inc., New York, 1982. 

6 N.Sh. Izmailian and Chin-Kun Hu, Physica A 254 (1998) 
198. 

7 T. Stosic, B.D. Stosic and LP. Fittipaldi, J. Mag. Mag. 
Mater. 177-181 (1998) 185. 

8 M.E. Fisher, in Critical Phenomena, Proceedings of the In- 
ternational School of Physics - Enrico Fermi, Course LI, 
edited by M.S. Green (Academic Press 1970). 

9 T.P. Eggarter, Phys. Rev. B 9 (1974) 2989. 



Finally, as T approaches zero, the multiplicative term 
(3 takes over and susceptibility diverges trivially as T _1 , 
for all generation levels. 

In conclusion, the exotic structure of the Cayley tree, 
with its infinite dimension and finite order of ramifica- 
tion, gives rise to rather unusual thermodynamic behav- 
ior. The overall scaling behavior of susceptibility is gov- 
erned by scaling of its amplitude, rather then distance 
from the critical point. In particular, cojitijary to the 
conclusions drawn from the previous worksBoci using the 
approximate formula (|l|), at Tc the susceptibility dis- 
plays an extremely weak singularity, with amplitude di- 
verging proportionally to the tree generation level n (as 
n increases to infinity in the thermodynamic limit), with 
the critical exponent 7 = 0. 
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